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The Cosmic Microwave Background Radiation (CMBR) holds infornration about almost all the 
fundamental cosmological parameters, and by performing a likelihood analysis of high precision 
CMBR fluctuation data, these parameters can be inferred. However, this analysis relies on as- 
sumptions about the initial power spectrum, which is usually taken to be a featureless power-law, 
P(k) oc k" B ~ . Many inflationary models predict power spectra with non-power law features. We 
discuss the possibility for detecting such features by describing the power spectrum as bins in k- 
space. This method for power spectrum reconstruction is demonstrated in practise by performing 
likelihood optimization on synthetic spectra, and the difficulties arising from reconstructing smooth 
f*"^ ■ features using discontinuous bins are discussed in detail. 
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O . PACS numbers: 98.70.Vc, 98.80.C, 98.80.-k 

o 

r\ ■ I. INTRODUCTION 

Fluctuations in the Cosmic Microwave Background Radiation (CMBR) were detected for the first time by the 
£f*s ', COBE satellite in 1992 [fLj. Subsequently it was realised that precision measurements of the fluctuation spectrum can 
J> ■ be used to infer almost all of the fundamental cosmological parameters [g]-f7| . The method most commonly used for 
determining parameters from the data is to maximise the likelihood function, £, over the space of model parameters 
to be determined, 6 = {f2, $7 m , Of, , Hq, t, . . .}. A large number of papers have dealt with with this issue in great detail. 
However, when calculating the theoretical CMBR power spectrum needed for the likelihood function, a necessary 
input is the initial power spectrum of fluctuations, usually assumed to have been produced during the inflationary 
epoch. 

There exists a plethora of different inflationary models, each having a specific prediction for the produced fluctuation 
spectrum. In the simplest slow-roll models, the spectrum is described by a simple featureless power-law P(k) oc fc" s_1 , 
where n s ~ 1 H. However, the power spectrum could easily behave differently, depending both on the specific form 
of the inflaton potential, V(<fi), and on the presence of different physical phenomena during inflation. A plethora of 
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different models predicting such behaviour exists M . Examples of such phenomena are the resonant production of 
+2 . particles during inflation, proposed by Chung et al. [Q, and the multiple inflation model discussed by Adams, Ross 
andSarkar 0. 

In almost all existing likelihood analyses the initial power spectrum is assumed to have the power-law form P(k) oc 
fc" s_1 , where n s is a constant [pH7|. This type of analysis restricts the parameter estimation in such a way that no 
non-power law features can be detected. Since there are so many different inflationary models, each with unique 
predictions, it is highly desirable to have a more model-independent way of estimating the initial power spectrum. 

Such a possibility has been discussed by Souradeep et al. jl^] and by Wang, Spergel and Strauss [[13],[L4|. In both 
treatments, the power spectrum was described as a featureless spectrum, but binned in fc-space. Using a Fisher 
matrix analysis, it was shown by Wang, Spergel and Strauss that it is possible to determine the power in each bin 
to reasonable precision, so that some features could be detectable. In the next section we discuss the prospects for 
detecting features in more detail, also using the Fisher matrix technique. Particularly, it is shown that there is a 
trade-off between the precision to which the power in each bin can be measured and the width of each bin. 

In Section III we go on to discuss in detail how such a likelihood analysis with a binned power spectrum is performed 
in practise. This analysis highlights some of the difficulties which can be expected. In general it is quite difficult to 
reconstruct a power spectrum with smooth features using a set of discontinuous bin amplitudes. If the binning is 
too coarse, the individual bins being comparable in size or broader than the power spectrum features, the spectrum 
reconstruction can easily give misleading results. For a reasonable reconstruction it is necessary to choose a binning 
which is significantly finer than the features to be detected. However, maximising the likelihood function is a highly 
non-linear optimization problem. If the number of free parameters in the fit, N = Ng + Nbins, is too large the 
computation time to retrieve the maximum likelihood becomes very long. Nevertheless, it is demonstrated that it is 
computationally feasible to use at least 20 bins in fc-spacc. 

Finally, Section IV contains a discussion of the results, with emphasis on how the techniques can be used on results 
from high precision satellite experiments, like MAP and PLANCK |lq| . 



II. FISHER MATRIX ANALYSIS 

CMBR temperature fluctuations are usually expressed in terms of spherical harmonics as 

AT 

— (M) = 5>Jm*Jm(M). (1) 
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From these a/ m coefficients it is possible to construct the power spectrum as 

Ci = (\a lm \ 2 ), (2) 

where the average in principle is an ensemble average. Such an ensemble average can obviously not be performed 
since we have access to only one realisation of the underlying distribution. However, in a universe which is isotropic 
the ensemble average can be replaced by an average over m- values for a given I JL7| . 

It is possible to estimate the precision with which the cosmological model parameters can be extracted from a given 
hypothetical data set. The starting point for any parameter extraction is the vector of data points, x. This can be 
in the form of the raw data, or in compressed form, either as the ai m coefficients or the power spectrum, Ci. Each 
data point has contributions from both signal and noise, x = jecmbr + ^noiso- If both signal and noise are Gaussian 
distributed it is possible to build a likelihood function from the measured data which has the following form Jig ] 

C{&)^e^>(- l -x^[C{Q)- l ]xY (3) 

where = (fi, ilf,, H , n,T, . . .) is a vector describing the given point in model parameter space and C(0) = (xx T ) is 
the data covariance matrix. In the following we shall always work with data in the form of a set of power spectrum 
coefficients, Ci. 

If the data points are uncorrelated so that the data covariance matrix is diagonal, the likelihood function can be 
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reduced to £ ex e~ x ' , where 



x 1 = ^ v ^' uuH ,,;r" y; , (4) 



2 V - ^ (,l>i,obs (w, theory) 



is a x 2 -statistics and iV max is the number of power spectrum data points |[L8| . 
The maximum likelihood is an unbiased estimator, which means that 

<9> = 9 . (5) 

Here 0o indicates the true parameter vector of the underlying cosmological model and (0) is the average estimate of 
parameters from maximising the likelihood function. 

The likelihood function should thus peak at ~ 0o, and we can expand it to second order around this value. The 
first order derivatives are zero, and the expression is thus 



x 2 







d6i 89 3 v * "ddidOj 



(0, - 6), (6) 



where i,j indicate elements in the parameter vector 0. The second term in the second derivative can be expected 
to be very small because (C; ]0 bs — Ci) is in essence just a random measurement error which should average out. The 
remaining term is usually referred to as the Fisher information matrix 

_ gy ^ ! dCidCi 

The Fisher matrix is closely related to the precision with which the parameters, 9i, can be determined. If all free 
parameters are to be determined from the data alone without any priors then it follows from the Cramer- Rao inequality 
HI that 

a(0i) = (F-% (8) 



for an optimal unbiased estimator, such as the maximum likelihood 20 
In general, the standard error on the Ci coefficients can be written as 
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Here, / s k y is the sky coverage and A oxp is a (Gaussian) experimental error. Notice that even for A oxp = 0, cr(Ci) =/= 0. 
This derives from the fact that the ensemble average in Eq. (|2|) has been replaced by an average over m-values, and 
is usually referred to as "Cosmic Variance" |TJJ . 

The CMBR is also predicted to be polarized, and this polarization power spectrum can in principle be measured. 
For scalar perturbations there are only three independent quantities: The temperature, C^t, the i^-ficld polarization, 
Ci^e, and the temperature-polarization cross-correlation, Ci t c- Then the Fisher matrix of Eq. ([?]) is instead [g(| 
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In this case the covariance matrix, Cov(Ci t x, C;.y), is a symmetric 3x3 matrix with the following elements 



Cov(Cut,Ci, t ) = 
Cov(C ltE ,Ci tE ) = 
Cav(jOi f c,C h c) = 
Cov(Ci,t,Ci,e) = 
Cov(C l>T ,Ci tC ) = 
Cov(C,E,C ltC ) = 



2 


(21 + l)/ sky 
2 


(21 + l)/ sky 
2 


(21 + l)/ sky 

2 


(21 + l)/ sky 
2 


(21 + l)/ sky 
2 



{Ci,t + A cxp! i)^ 



(CS 



l.E 



icxp,2j 



v-'i.e + (C/,t + A exP] i)(Cj i _E + A exP]2 )) 



r 2 



(11) 



(21 + l)/ sky 



Cl,c(Cl,T + A eX p,l) 

Ci,c(Ci,e + A exPj2 ), 



where A oxp yi are again experimental experimental errors related to pixel noise and beamwidth pOfl , 

For the purposes of the present paper we assume that / sky = 1, A cxp .i = and Z max = 1500, corresponding to 
a full-sky survey up to l max = 1500, limited only by cosmic variance. If the temperature power spectrum alone is 
considered this is not too far from what can be expected with MAP, whereas PLANCK will likely measure both 
temperature and polarization to this accuracy. 

As the free cosmological parameters we use the matter density, r2 m , the cosmological constant, fl\, the baryon 
density, fif,, the Hubble parameter, H$, the optical depth to reionization, r, and the overall normalization, Q. Instead 
of using the spectral index of a featureless power-law spectrum, n s , as the last free parameter, we bin the spectral 
index in fc-space. In practise we use N bins of equal size in log(fc), from k = 1.27 x 10 -5 Mpc -1 to k = 0.25 Mpc -1 
(this covers the entire range of /c-space which is visible in the CMBR) . In each bin the power spectrum is assumed to 
follow a power-law, P(k)i oc k n "^~ 1 . Thus, the parameter vector is 



& = {n m ,n A ,n b ,H ,T,Q,n s (l),...,n 8 (N)} 



(12) 



As the reference model around which to calculate the Fisher matrix, we take the standard CDM model with, fl = 
fl m = 1, Oft = 0.05, Hq = 50 kms _1 Mpc , and r = 0. The reference model has a power-law initial spectrum with 
n s (i) = 1 for all i. 





2.0 




1.5 




1.0 




0.5 




0.0 




2.0 




1.5 




10 



10 " 10 
k (Mpc" 1 ) 



10 



FIG. 1. The precision with which the spectral indices in each bin can be measured, shown for three different numbers of bins 
(TV = 1, 5, 20), for the case where only the temperature anisotropy can be measured. Shown are the ler error bars around the 
central value which in all cases is n a = 1. 



In Fig. 1 the precision with which n s (i) can be measured is shown for different numbers of bins, for the case where 
only the temperature power spectrum can be measured. Not surprisingly, the precision with which the effective spectral 
index in each bin can be measured depends strongly on the number of bins. For only one bin, An s /n s — 8.21 x 10 , 
whereas for 20 bins, the smallest An s /n s is 3.49 x 10~ 2 . Thus, there is a trade-off between the resolution in fc-space and 
the resolution in n s . Detecting a low-amplitude narrow feature is not possible, whereas either a broad low-amplitude 
or a narrow high-amplitude feature should be detectable. In agreement with Wang, Spergel and Strauss [ [J3|]l4l ] we 
find that the CMBR data is most sensitive in the range k ~ 10~ 3 — 10 _1 Mpc . If polarization can also be measured, 
the precision with which the n s (i) can be measured increases. Fig. 2 shows the same as Fig. 1, but with the inclusion 
of polarization. Indeed, the precision on the individual bin amplitudes increases by a large factor. For only one bin, 
An s /n s = 2.90 x 10~ 3 , whereas for 20 bins, the smallest An s /n s is now 9.49 x 10~ 3 . 
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FIG. 2. The precision with which the spectral indices in each bin can be measured, shown for three different numbers of bins 
(TV = 1,5,20), for the case where both temperature and polarization can be measured. Shown are the la error bars around 
the central value which in all cases is n s — 1. 



Next, a fundamental question is how much the ability to determine the other fundamental parameters depends on 
the number of bins in fc-space. Figs. 3 and 4 show the precision in measuring the other cosmological parameters as a 
function of N, the number of bins. Fig. 2 is for the case where only the temperature anisotropy can be measured. For 
all the cosmological parameters, except r, there little degeneracy between these parameters and the power spectrum 
indices, meaning that there is little loss of ability to pin down these parameters. From N = 10, the resolution begins 
to decrease, although quite slowly, and at N = 30 it is for instance down by about 40% for £V Clearly, the spectral 
indices, n s (i), are not very degenerate with the other cosmological parameters. The reason is that the discontinuous 
binning in the power spectrum is very hard to mimic by continuous changes in other parameters. This finding is in 
agreement with Ref. |13| , |l4|, where such a binning was also used. However, there is a marked difference compared 
with the findings of Ref. ]12| . In this work, smoothness of the power spectrum was enforced by using smooth shape 
functions peaked at a given point in fc-space instead of discontinuous binning. These authors find that the uncertainty 
in the other parameters increases very rapidly with N, rising by more than a factor of 10 at JV = 30 E3|. However, 
using smooth functions to describe the power spectrum could indeed be expected to be able to mimic changes in 
other parameters to a much higher degree than our approach. So it is not surprising that a much higher degree of 
degeneracy is found. On the other hand, as will be seen in the next section, using a discontinuous binning has the 
very clear disadvantage that it can be difficult to achieve reasonable likelihood fits if the number of bins is too small. 



CD 

b 



0.10 
0.08 

0.06 
0.04 h 
0.02 



0.00 



CT (n A ) 

a(n b )/0.05 

h /0.5 

t/5 

/ 



../ 



/ 



No. /c-bins 



FIG. 3. The precision with which the other cosmological parameters, fl m ,^A,iib and Ho, measured as a function of the 
number of bins in fc-space, for the case where only the temperature anisotropy can be measured. The precision is quantified in 
terms of the standard error, a. 

In line with this discussion, it was shown by Kinney [|l5| that smooth power spectrum features can mimic variations 
in the cosmological parameters almost exactly (to within cosmic variance). In this work 75 bins in fc-space were 
used and the power spectrum smoothcncd by cubic spline interpolation, leading to a very large degree of degeneracy 
between the fc-space power spectrum and the other cosmological parameters. 

The one exception to the above is the optical depth to reionization, r. Here there is a very large degree of degeneracy, 
which increases rapidly with the number of bins. This effect was also described by Wang, Spergel and Strauss [13[. 
It happens because the effect of reionization is to suppress the amplitude at small scales by a factor e _2r . This can 
be quite easily mimicked by suppressing the fc-space power spectrum at small scales. 

However, Fig. 4 shows the case where polarization can also be measured. In this case, the degeneracy between r 
and the fc-space spectrum is broken, and the ability to determine the value of t is increased by a large factor. For the 
other cosmological parameters, the precision is also increased, but not as dramatically. 
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FIG. 4. The precision with which the other cosmological parameters, Q m ,QA,Qb and Ho, measured as a function of the 
number of bins in fc-space, for the case where both temperature and polarization can be measured. The precision is quantified 
in terms of the standard error, a. 



III. NUMERICAL SPECTRUM RECONSTRUCTION 



Having done this initial estimate of how precisely the initial power spectrum can be estimated, it is important to 
see how such a power spectrum reconstruction works in practise. In order to investigate this, we have produced 10 
synthetic spectra based on a underlying theoretical model. These spectra are calculated assuming Gaussian errors 
given by cosmic variance only. In order to simulate a possible feature in the power spectrum we introduce a Gaussian 
"bump" , so the spectral index has the form 



(log(fc*) - log(fc ,*)) : 



n s = l-Acxp{- ^^*> "&™*» \ (13) 

where fe* = fc/lMpc~ . This bump is characterised by three parameters: A, the amplitude, fco,*, the position in 
fc-space, and a, the width. For the underlying cosmological model we have chosen the same standard CDM model as 
in the last section, so the model parameters are: fl — fl m = 1, Qj, = 0.05, Ho = 50 kms -1 Mpc -1 , A — 0.7, a = 0.1 
and logfco,* = — 2- Thus, the bump has been placed closed to the region where the CMBR data is most sensitive. 
Note that in this section we use only information related to the temperature power spectrum, not the polarisation 
spectrum. 

On each of these synthetic spectra we try reconstructing the power spectrum using the logarithmic binning method 
introduced in the previous section. The likelihood optimization algorithm is based on the simulated annealing principle 
described in Ref. J21fl . This method has the advantage that it is very fast for highly non-linear optimization over 
many-dimensional parameter spaces, exactly the problem that is faced here. In Fig. 3 we show the values extracted 
by the parameter reconstruction for three different numbers of bins, N — 5, 10 and 20. 

In the N = 5 case the individual bins are broader than the bump to be reconstructed. This shows up in the fact 
that no good fit is achieved ((Xmi n )/d.o.f = 12.7), and the spectrum reconstruction is quite poor. The reconstructed 
spectrum does show some evidence of a lower spectral index around the position of the bump, but the magnitude is 
not correct. For N = 10 the reconstruction is already much better and reproduces the actual slope of the underlying 
spectrum. However, (x^i^/d.o.i = 8.06, so the fit is not very good in this case either. In the last case where N = 20 
the reconstructed spectrum again follows the underlying one reasonably well. However, at N = 20 the uncertainty on 
the individual bin amplitudes is already significantly larger than for N — 10. By going to an even higher number of 
individual bins, the spectrum bump will be smeared out beyond recognition. The x 2_n t is significantly better than 
for N = 10 ((x, 2 n ; n )/d.o.f = 2.84), although far from what should be expected from a "good" fit ((x, 2 n i n )/d.o.f ~ 1). 
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FIG. 5. The reconstructed power spectrum from likelihood maximisation over 10 synthetic spectra. The boxes show the lcr 
error intervals. The full line shows the underlying fc-space power spectrum used to generate the synthetic spectra. 

The reason for these very high x 2 values is that it is quite hard to mimic a smooth feature in the power spectrum 
with a sequence of discontinuous bins. Thus, even though it is indeed possible to map out the shape of the power 
spectrum, it is not possible to achieve a good fit in terms of the likelihood without letting N — ► d.o.f . However, going 
much beyond N — 20 increases the uncertainty in the individual bin amplitudes and doing the likelihood maximization 
is already quite demanding numerically at N = 20. 

Finally, with a large number of bins (10 or 20) the reconstruction at scales smaller than the bump (large k) is 
consistently predicting a too low spectral index. This effect is worrying and to investigate whether it is a generic 
feature of any discontinuous binning method, or specific to our choice, where each bin has a "tilt", P(k)i oc fc™" 1 ' -1 , 
we also do the reconstruction using the same method as Wang, Spergel and Strauss p3[ . Here, the power in each bin 
is assumed to have constant amplitude over the entire bin so that P{k)i = Aj. We show the reconstruction for the 
case of N = 20 in Fig. 6. For this case, the reconstruction of the power spectrum bump is actually better than for 
the "tilt" spectrum reconstruction. There is no trace of the spuriously low spectral index found at high k. However, 
the x 2 is substantially worse for this method, with (x^ lin )/d.o.f = 3.45. From this, a robust conclusion is that for 
discontinuous binning the obtainable fits are quite poor in terms of x 2 , regardless of the binning method. However, it 
is possible to detect features in the spectrum by this method. At least the tilted binning method can show spurious 
effects, so one should be careful about making strong conclusions without going through several different methods of 
spectrum reconstruction. 
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FIG. 6. The reconstructed power spectrum from likelihood maximisation over 10 synthetic spectra, for the case of N = 20. 
The boxes show the lcr error intervals. The full line shows the underlying fc-space power spectrum used to generate the synthetic 
spectra. 



Just for comparison we have also done a spectrum reconstruction where the free parameters are those describing 
a smooth scale invariant power-law spectrum (n s = 1) with a Gaussian bump: A, &o,* and a, instead of using n s (i). 
The values extracted from likelihood maximisation on the 10 synthetic spectra are given in Table I. Clearly, in this 
case, the true underlying spectrum is recovered to very good accuracy because the same functional form is used in 
the fit as in the underlying power spectrum. Also the x 2 values are completely consistent with what is expected from 
a good fit. 

So using smooth fitting functions could have the advantage of being better able to fit smooth features in the 
underlying initial power spectrum. However, the discontinuous binning method has the advantage of being very 
model independent. So a possible way to proceed would be to initially try the discontinuous binning method, and 
only for refinements go to more complicated fitting functions. 



IV. DISCUSSION 

We have discussed the possibility for measuring the shape of the primordial power spectrum by using precision 
measurements of CMBR fluctuations. The method employed was to parametrise the spectrum as a power-law, binned 
in fc-space. An analysis of the obtainable precision in measuring the power spectrum, using a Fisher matrix analysis, 
showed that it should be possible to detect features that are sufficiently broad by using CMBR observations. 

However, there are some practical problems involved in this spectrum reconstruction. If the power spectrum is 
parametrised by binning it in fc-space, it is very difficult to mimic any smooth features in the underlying power 
spectrum. Especially if the feature is comparable in width to the individual bins the reconstructed spectrum can show 
spurious features. Nevertheless, by going to a sufficiently fine binning it is possible to map out the general shape of 
the power spectrum, although still very difficult to get a good fit in terms of the likelihood. 

The fact that it is possible to reconstruct the underlying fc-space power spectrum from CMBR observations opens 
up the possibility of probing physics at the time of spectrum formation by observing the universe at CMBR formation 
(T ~ 4000K). As was for instance discussed in Refs. p3Lp2[, if the primordial fluctuations are produced during the 
inflationary epoch, it is possible to reconstruct the inflationary potential, V(<j)), from the fc-space spectrum. 

However, in addition to that it will be possible to detect non-power-law features in the power spectrum using this 
discrete binning method. Such features are for instance predicted to occur in the multiple inflation model of Adams, 
Ross and Sarkar pT| . Here, the effect of symmetry breaking during the inflationary period was discussed. If such 
spontaneous symmetry breaking occurs along flat directions, short inflationary periods will result, leaving distinct 
non power-law features in the power spectrum. Another possible mechanism for producing features is the resonant 
production of particles during inflation discussed by Chung et al. pi . 

Note that the CMBR is mainly sensitive to power spectrum at fc ~ 10 -3 — 10 _1 Mpc^ 1 . Any features outside this 
region in fc-space would not be detectable, so the CMBR can at most probe a very small region of the inflationary 
period. However, as discussed for instance in Ref. JRjjlJ], by also using data from large scale structure (LSS) 
surveys such as the Sloan Digital Sky Survey it is possible to increase the region of fc-space that can be probed. 
LSS surveys are sensitive to smaller scales than the CMBR (for instance the Sloan survey should probe the region 
fc ~ 10~ 2 — 1 Mpc -1 ), and by combining MAP/PLANCK with such surveys it should be possible to have sensitivity 
in the region fc ~ 10~ 3 — 1 Mpc -1 . In addition to this, including data from large scale surveys can break some of the 
cosmological parameter degeneracies H. For instance, using CMBR data alone it is impossible to determine either fl m 
or 51a separately with high precision. The parameter that can be measured is the total energy density £1 = Q m + Q\. 
However, the LSS surveys are sensitive to a different combination of £l m and SIa, and by adding this information the 
degeneracy between them can be broken H] . 
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TABLE I. The recovered parameters from the fit to a 
Gaussian bump overlayed on a scale invariant spectrum. 10 
synthetic spectra were used for the reconstruction. 



Parameter 



found 



Amin 
Q 

logfco, 
A 



expected 



1516.7 ±50.4 


1496 ± 54.7 


0.1004 ±0.0035 


0.1 


-2.0004 ± 0.0059 


-2 


0.6931 ±0.003 


0.7 
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